Computational Ghost Imaging 
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Ghost- imaging experiments correlate the outputs from two photodetectors: a high spatial- 
resolution (scanning pinhole or CCD camera) detector that measures a field which has not interacted 
with the object to be imaged, and a bucket (single-pixel) detector that collects a field that has in- 
teracted with the object. We describe a computational ghost-imaging arrangement that uses only 
a single-pixel detector. This configuration affords background-free imagery in the narrowband limit 
and a 3D sectioning capability. It clearly indicates the classical nature of ghost-image formation. 

PACS numbers: 42.30.Va, 42.50.Ar, 42.25.Kb 



Ghost imaging is the acquisition of object information 
by means of photocurrent correlation measurements. Its 
first demonstration utilized a biphoton source, thus the 
image was interpreted as a quantum phenomenon owing 
to the entanglement of the sourcephotons [1]. Exper- 
imental [2, 3] and theoretical 0, [5| work later demon- 
strated that ghost imaging could be performed with pseu- 
dothermal light. Whereas the biphoton requires a quan- 
tum description for its photodetection statistics, pseu- 
dothermal light can be regarded as a classical electro- 
magnetic wave whose photodetection statistics can be 
treated via the semiclassical (shot-noise) theory of pho- 
todetection [6]. This disparity has sparked interest in 
the physics of ghost imaging BSSIl^. Recently [Tl|, 
we established a Gaussian-state analysis of ghost imag- 
ing that unified prior work on biphoton and pseudother- 
mal sources. Our analysis indicated that ghost-image 
formation is intrinsically due to classical coherence prop- 
agation, with the principal advantage afforded by the 
biphoton state being high-contrast imagery in the wide- 
band limit. Other recent work [12, 13], however, has 
claimed that pseudothermal-light ghost imaging is a fun- 
damentally quantum phenomenon — one of nonlocal two- 
photon interference — that cannot be explained in terms 
of intensity-fluctuation correlations. In this Letter we 
shall show that ghost imaging can be accomplished with 
only one detector, viz., the bucket detector that collects 
a single pixel of light which has been transmitted through 
or reflected from the object. As only one light beam and 
one photodetector are required, this imaging configura- 
tion cannot depend on nonlocal two-photon interference. 
Moreover, it affords background- free imagery in the nar- 
rowband limit and a 3D sectioning capability. 

We begin with a review of pseudothermal-light lens- 
less ghost imaging, based on 11 , using classical elec- 
tromagnetism and semiclassical photodetection theory. 
Consider the setup shown in Fig. 1. Here, Es{p,t)e~^^^^ 
and Eji{p^ t)e~*'^°^ are scalar, positive frequency, classical 
signal (S) and reference (R) fields that are normalized to 
photon- units and have center frequency ujq- They are z- 
propagating with p being the transverse coordinate with 
respect to their optical axes. More importantly, they 
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FIG. 1: (Golor online) Pseudothermal ghost-imaging setup. 



are the outputs from 50/50 beam splitting of E{p^t), a 
continuous-wave (cw) laser beam that has been trans- 
mitted through a rotating ground-glass diffuser. The sig- 
nal and reference undergo quasimonochromatic parax- 
ial diffraction over L-m-long free-space paths, yielding 
measurement-plane fields [14] 



dp^E^{p\t-L/c)^ — , (1) 



where (i^m) = (1,5*) or (2,i?), c is the speed of light, 
and ko = u;o/c = 27r/Ao. The field Ei{p^t) illuminates a 
shot-noise limited pinhole photodetector centered at Pi 
with sensitive region p e Ai. The field E2{p,t) illu- 
minates an amplitude-transmission mask T{p), located 
immediately in front of a shot-noise limited bucket pho- 
todetector with sensitive region p e The product 
of the photocurrents from these detectors is time aver- 
aged to estimate their ensemble-average cross correlation, 
C{pi). This process is repeated, as Pi is scanned over the 
plane, to obtain the ghost image of the object's intensity 
transmission |T(p)p [lH. We then have that 

C(pi) = ^Jdriqr]iPi{t-ri)hi{ri) 

X J dT2 qV2P2{t - T2)h2{T2)^ . (2) 
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In this expression: q is the electron charge; rj£ is the 
quantum efficiency of the pinhole = 1) and bucket 
{£ = 2) detectors; Pi{t) = j^^dp\Ei{p,t)\^ and = 
jj^^dp\E2{p^t)\^\T{p)\^ are the photon fluxes imping- 
ing on detectors 1 and 2; and hi{t) is a baseband im- 
pulse response representing the finite response time of 
the photodetector i. We shall assume that the pinhole 
is small enough that Pi(t) ^ Ai\E{pi^t)\'^ , where Ai 
is that detector's photosensitive area. In terms of the 
intensities (photon-flux densities) t) = |£^i(p, 

and hip^t) = \E2{p,t)\'^\T{p)\'^ that illuminate the pho- 
todetectors we then have that 

{Pi{h)P2{t2)) = Ar{h{p,,h)) I dp{l2{pM)) 

^ Ai [ dp(A/i(pi,ti)A/2(p,t2)), (3) 

where A/£(p, t) = Ii{p^ t) — (/^(p, t)) is the intensity fluc- 
tuation. The flrst term on the right gives rise to a fea- 
tureless background, while the second term leads to the 
ghost image, as we now show for a Gaussian-Schell model 
of pseudothermal illumination. 

Let E[p^t) be a zero-mean, cross-spectrally pure [l6| . 
complex- valued Gaussian random process that is charac- 
terized by its phase-insensitive correlation function 

{E%p^M)E{p2M)) = K{p,,p^)R{t2 - ti), (4) 

where R{0) = 1. Using Em{p,t) = E{p,t)/V2, for m = 
5', R, and Eq. ^ we have that {Ei{p, t) : £ = 1, 2}, is a 
pair of completely correlated, zero-mean, complex-valued 
Gaussian random processes that are characterized by the 
following auto- and cross-correlation functions: 

(^;(Pl,ti)^,(p2,fe)) = (^r(Pi,ti)^2(P2,fe)) (5) 



with aL = 2L/kopo and pL = 2L//co<^o- Physically, 
and pz are the intensity radii and coherence radii of the 
flelds at 2: = and 2; = L, so that the preceding behav- 
ior represents the familiar situation for partially coherent 
light in which the far-fleld intensity radius is controlled 
by the source's coherence radius and the far-field coher- 
ence radius is controlled by the source's intensity radius. 

Suppose that the photodetector impulse responses 
h£{t) have response times that are much shorter than the 
field's coherence time. To, i.e., we are in the narrowband 
regime. The Gaussian-Schell model source then leads to 



X / dpe-'l^l'/'^'^-l^i-^l'/^'H^lp)!', (10) 
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where we have assumed that J dthi{t) = 1, and 



Co(pi) = q r]imAi[ — 2 ' ^ 
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When T(p) and pi are space limited to a radius much 
smaller than about the origin, C{pi) is comprised of 
a constant background term. 



Co = q ViV2( — 2 



dp\T{p)f 



(12) 
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plus the ghost-image term 



K\p^,p^)R{t2-h), (6) Ci(pi) =^'7/1 7/2^1 



for £ = 1,2. Given an explicit K{pi,p2) for the spa- 
tial auto-correlation function of E{p^t), it is a simple 
matter, in principle, to calculate K'{p^^p2) via standard 
coherence-propagation theory 0. We then have that 
{h{p,t)) = K'{p,p) and {h{p,t)) = K'{p,p)\T{p)\\ 
More importantly, the moment-factoring theorem for 
Gaussian random processes [17] implies that 

(A/i(pi,ii)A/2(P2,i2)) = \K'{p^,p^)\^\R{t2-ti)\^ 

X \np,)\\ (7) 

In the far field (when koaopo/2L <C 1) the Gaussian- 
Schell model correlation function for E{p^t)^ 

2P 



K{pi,P2) = 
with po <C ao, yields 



-(|pir+|P2r)/«o-iPi-P2r/2p^ 



(8) 
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(9) 



/ dpe-l^i-^l'/^' |T(p)|^ 

JA2 

(13) 

Equations (p!2|) and (p!3|) summarize the key elements of 
pseudothermal ghost imaging. Within an object-plane 
region whose spatial extent is small compared to XoL/po, 
we obtain a pseudothermal-light ghost image with spatial 
resolution ^XoL/ao [l^ that is embedded in a featureless 
background [19] . The background term can be eliminated 
by employing a zero-frequency (DC) block between one 
or both of the photodetectors and the correlator shown 
in Fig. 1, as done in the experiment reported in [lot . 

With the preceding analysis in hand, it becomes a sim- 
ple matter to walk our way through to a single-pixel 
ghost imager. First, rather than use laser light trans- 
mitted through a rotating ground glass as the source of a 
narrowband, spatially- incoherent £^(p, t), let us employ 
the configuration shown in Fig. 2. Here, we transmit a 
cw laser beam through a spatial light modulator (SLM) 
whose inputs are chosen to create the desired coherence 
behavior. Specifically, we assume an idealized SLM con- 
sisting of dx d pixels arranged in a (2M + 1) x (2M + 1) 
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FIG. 2: (Color online) SLM ghost- imaging setup. 

array with 100% fill factor within di D x D opaque pupil, 
where D = (2M + l)d and M > 1. We use this SLM to 
impose a phase (j)nm {t) on the light transmitted through 
pixel (n.m)^ with {e*'^'^'^^^) : —M < n^m < M} be- 
ing independent identically-distributed (iid) random pro- 
cesses obeying (e^^--^^)) = and (e^([^--(^2)-0ife(^i)]^ = 
^jn^/cme~'*^~*^'/^°, where the coherence time Tq is long 
compared to the response times of the h£{t) (2q| . 

In the far field, i.e., when k^dD / L <C 1, the preced- 
ing E{p,t) leads to Ei{p^t) and E2{p^t) that are zero- 
mean fields whose correlation functions satisfy Eq. (|6]) 
with R{r) = e-l^l/^o and 



= f (ill 
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(14) 



Although it is not a Gaussian-Schell form, the preced- 
ing spatial correlation function has an intensity width 
r^XoL/d dLiid a coherence length ^XqL/D^ behavior which 
is qualitatively similar to what we saw earlier if we iden- 
tify d ^ po and D ^ qq. Furthermore, Central Limit 
Theorem considerations imply that ^i(p, t) and ^2(p, t) 
may be taken to be jointly Gaussian. Hence our Fig. 2 
configuration will produce a ghost image of spatial res- 
olution XqL/D within a spatial region of width XoL/d 
embedded in a background by virtue of 

= q%r]2ArK\p,,p,) f dp K\p, p)\T{p)f 
+ q^mmAi I dp\K\p^,p)\^\T{p)\\ (15) 

As before, the background term can be suppressed, if 
desired, by means of a DC block. 

To realize the Fig. 2 ghost imager we could use noise 
generators to drive the SLM in a way that approxi- 



mates the preceding statistics. It is more interesting, 
for what will follow, to suppose that deterministic driv- 
ing functions are employed to achieve the same objec- 
tive. What we want at the SLM's output is a narrow- 
band, zero-mean field whose spatial coherence — inferred 
now from a time average, rather than an ensemble av- 
erage, because there is no true randomness — is limited 
to field points within a single pixel. Sinusoidal modula- 
tion, (j)nm{t) = ^ cos(r^o + A^nm^), with different /^Vtnm 
for each pixel will work. Let (•)t^ denote time averaging 
over the T^-sec interval employed in obtaining the ghost 
image. We have that (e*'^^"^'^^^)T^ Jo(<^>) ~ 0, where Jo 
is the zeroth-order Bessel function of the first kind, when 
{VLq + /\Vtnm)Ta > 27r and <l> > tt. With the additional 
condition \/^Vtnm\ < ^0 we have {e^^'t>r^rr^{t)-c|,jk{t)]^^^^ 
Jo(2<l>) ^ 0, unless j = n and k = m. Furthermore, the 
output field will satisfy our narrowband requirement if 
the modulation periods 27r/{Qo + Antrim) are all much 
longer that the response times of the h£{t). Thus, this 
deterministically-modulated source will also yield a ghost 
image of spatial resolution XqL/D within a spatial region 
of width XoL/d embedded in a background that can be 
suppressed by means of a DC block. 

At this point, the notion of computational ghost 
imaging — in which we only use the bucket detector — is 
easily understood, see Fig. 3. We use deterministic mod- 
ulation of a cw laser beam to create the field E2{p^t) 
that illuminates the object transparency, and, as usual, 
we collect the light that is transmitted through the trans- 
parency with a bucket (single-pixel) detector [21]. Know- 
ing the deterministic modulation applied to the original 
cw laser beam allows us to use diffraction theory to com- 
pute the intensity pattern, /i(pi, t), that would have illu- 
minated the pinhole detector in the usual lensless ghost 
imaging configuration. Furthermore, we can subtract the 
time average of this intensity, in our computation, and 
obtain the equivalent of AIi{pi^t). To distinguish these 
computed values from actual light-field quantities, we will 
denote them Ii{pi,t) and AIi{pi,t)^ respectively. The 
time average correlation function. 



A(7(pi) 



dTi g7?i Ai A/i (pi , t - Ti ) (ri ) 



j dT2 qr]2P2{t - T2)h2{T2) 



(16) 



will then be a background-free ghost image — with spatial 
resolution XqL/D over a spatial extent of width X^L/d — 
akin to what would be obtained from pseudothermal 
ghost imaging with d ^ po^ D ^ ao^ and a DC block 
applied to the pinhole detector. Now, because only one 
photodetector has been employed, it is impossible to in- 
terpret this computational ghost image as arising from 
nonlocal two-photon interference. 

In addition to obviating the need for a high spatial- 
resolution detector in ghost-image formation, at the ex- 
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FIG. 3: (Color online) Computational ghost- imaging setup. 

pense of the computational burden associated with the 
free-space propagation calculation for the field £'(p, t), 
computational ghost imaging permits 3D sectioning to 
be performed. To see that this is so, we return to the 
Fig. 1 ghost-imaging configuration with pseudothermal 
light and inquire about its depth of focus,. In other 
words, how badly is the ghost image blurred if the object 
is at z = L but the pinhole detector is at z = L + AL? 
This question is easily answered. Equation ([7j) becomes 

(A/i(pi,ti)A/2(p2,i2)) = 

\K"{p„p,)\'\R{t2-h-AL/c)\^\T{p,)\^ (17) 



where 



K"{p,,p2)= dpK'ip,p2 



27rAi ■ 
(18) 

As a result, the ghost-image term from Eq. ([13]) 
takes the same form with pL replaced by p'^ = 
pL\/^ + (^^/^opi)^5 so that the focal region is |AL| < 
^oPl = 4L^//cotto- near-field of the pre-diffuser 

laser beam, i.e., when koal/^L 1, the focal region is 
a very small fraction of the source-to-object path, be- 
cause \AL\/L = AL/kottQ <C 1 as reported for the ex- 
periments in [12]. This limited depth of focus has the 
following implications when a range-spread opaque ob- 
ject is imaged in reflection. The pseudothermal ghost 
imager can only image one focal region at a time. How- 
ever, because the computational ghost imager can pre- 
compute A/i(p^,t) for a wide range of propagation dis- 
tances, the same bucket-detector photocurrent can be 
correlated with many such AIi{pi, t) to perform 3D sec- 
tioning of the object's reflectance. Of course, this sec- 
tioning further increases the computational burden, but 
this burden can be handled off-line, and, for a given SLM 
and its associated modulation waveforms, the same pre- 
computed AIi{pi^t) can be used for all ghost images 
formed using that system. 

In conclusion, we have shown that ghost imaging can 
be performed with only a bucket (single-pixel) detector 
by precomputing the intensity fluctuation pattern that 
would have been seen by the high spatial-resolution de- 
tector in lensless ghost imaging. This computational 
ghost imager yields background-free images whose res- 
olution and field of view can be controlled by choice of 
spatial light modulator parameters, and it can be used to 
perform 3D sectioning. Finally, the computational ghost 



imager underscores the classical nature of ghost-image 
formation. In particular, as only one light beam and one 
detector are employed, it is impossible to consider the 
computational ghost image as arising from nonlocal two- 
photon interference, as has been argued in [13, [13] for 
the pseudothermal case. 
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